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Abstract 

We study with lattice Monte Carlo simulations the relation of global 0(2) symmetry 
breaking in three dimensions to the properties of a geometrically defined vortex loop 
network. We find that different definitions of constructing a network lead to different 
results even in the thermodynamic limit, and that with typical definitions the perco- 
lation transition does not coincide with the thermodynamic phase transition. These 
results show that geometrically defined percolation observables need not display uni- 
versal properties related to the critical behaviour of the system, and do not in general 
survive in the field theory limit. 
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1 Introduction 



In a classic work [|TJ, Banks, Kogut and Myerson showed that the partition function of 
a particular three-dimensional (3d) spin model with global 0(2) symmetry, called the 
Villain model, can equivalently be represented as a partition function of a dual theory 
in which the spin configurations are integer-valued and sourceless and can therefore 
be represented as configurations of closed loops. Moreover, they suggested that at the 
transition point where the global 0(2) symmetry gets restored, infinitely long loops 
become important. The loops of the dual theory can be identified with vortex lines 
of the original theory, and therefore this result suggests that also in other theories 
that contain topological defects, they might play an important role in determining the 
properties of the phase transition. As a consequence, many problems in condensed 
matter physics and cosmology could be simplified by concentrating only on the vortex 
loop degrees of freedom (see, e.g., and references therein). 

In many realistic cases, this intuitive picture can be realized in the precise sense that 
the vortex tension T defined as the energy per unit length of an isolated vortex vanishes 
at the transition point 0, [|, |5|. Near the transition point, vortex loops behave as world 
lines of particles in 2+ld, and following the Villain case || one can then conjecture 
that they could be described with an effective scalar field of mass T ||, |. However, 
it is rather difficult to measure the tension T in practice || ||, |J, which limits the 
practical uses of this approach. 

Therefore, one may seek for alternative ways of studying the vortex degrees of free- 
dom. A seemingly natural approach is to decompose every spin configuration generated 
in a lattice Monte Carlo simulation into a number of closed vortex loops. The hope is 
then that the transition could be identified with a non-zero probability of finding vor- 
tex loops that extend through the whole system, a phenomenon which is often called 
percolation. 

Percolation has been used to study phase transitions in various different theories. 
Examples of this in condensed matter physics include the 3d XY (or Heisenberg) 
model |H| [11], |T2"|| , and in cosmology the corresponding global 0(2) field theory [p~3|j , 
both assumed to represent some features of the gauged Ginzburg-Landau theory where 
the loops would be interpreted as Nielsen-Olesen vortices. In particle physics, a similar 
approach has been used to give a physical interpretation for the confinement phase 
transition in non-Abelian gauge theories |TJ| and to define an order parameter in the 
crossover regime of the electroweak theory |TJ[ . 



Given the large range of physical applications for the percolation picture, it is im- 
portant to check whether it really works. Experience from the Ising model where a 
different kind of percolation may occur, related to clusters instead of vortex loops, 
has shown |16| that one has to be quite careful with this interpretation (for a recent 
discussion, see fl7fl ). In the present paper we study with high precision lattice Monte 
Carlo simulations the geometric properties of a properly defined vortex loop network 
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near the thermodynamic phase transition in the XY model and show that in all of 
the cases considered, the percolation point fails to coincide with the thermodynamic 
singularity. Nevertheless, a description in terms of vortex degrees of freedom can still 
be useful in many other contexts, if used with caution. 



2 Model 

The 3d XY model is defined by 

Zxy = rn x d0 x exp(/?££ cos (0*+i - X )). (2.1) 

J -t x i=l 

It belongs to the same universality class as the Villain model where the original con- 
sideration involving vortex loop degrees of freedom was carried out Jl|. However, it is 
calculationally simpler, its critical properties are known very accurately, and it is more 
directly related to continuum field theories such as the Ginzburg-Landau model. 

Let us recall that the only parameter point of the XY model corresponding to a 
continuum field theory is the 2nd order transition point at /3 — j3 c . Thus for obtaining 
results free from lattice artifacts, one should focus on the vicinity of f3 c . From the 
statistical mechanics point of view, this is the region where physical observables show 
universal critical behaviour. 



3 Method 

In order to analyse accurately the universal behaviour close to (3 C , we shall employ the 
finite size scaling method [Tj|. We recall that the idea of the method is to first locate 



the infinite- volume critical point of the system. One then measures finite- volume values 
for various observables at that point, and determines from them the critical exponents. 
We will here address particularly two questions: 

1. Does the infinite- volume extrapolation of the percolation critical point coincide 
with the thermodynamic critical point? 

2. Sitting at the thermodynamic critical point, do percolation related observables 
show critical behaviour? If so, do the critical exponents match any of the known 
thermodynamic ones? 

As a benchmark with which to compare, let us recall the known thermodynamic 
properties of scalar models with global 0(2) symmetry. In the XY case, the infinite 



volume critical point is at []19 



(3 C = 0.454165(4). (3.1) 



2 



Let us denote 

E = -£]Tcos(# x+i -0 X ), M = £cos(0 x ), (3.2) 

x i=l X 

where we have assumed the symmetry breaking direction to be 9 = 0. For any operator 
Oi, one can define the related susceptibility 

C i = iV 3 ((O i -(O i )) 2 ), (3.3) 

where N 3 is the volume of the system. For energy (E) and magnetization (M) like 
observables at /3 C , the susceptibilities diverge as 

C E oc N a ' u , C M oc N^ v . (3.4) 

Here v is the correlation length critical exponent. Magnetization itself scales at the 
critical point according to 

(|M|) ocN-V". (3.5) 

Finally, corrections to scaling are determined by an universal exponent u; for instance, 
the finite volume apparent critical point determined from C is assumed to behave as 

^-^ = ^ + 0,^ + .... (3.6) 



The known results for the exponents appearing here are [19, 20, 21 



a (3 7 v lo a/v /3/v 7/V 
035 L32 67 79 ^002 52 TW 



(3.7) 



4 Observables 

In order to locate the vortices in a given field configuration, we first define 

= [0x+i - 0*]*, (4.1) 

where [...]„ means (...) mod 2ir with a result in the interval (— 7r, 7r]. The winding 
number of a plaquette is then defined as 



n 



2tt 



(^,x + ^},x+i — ^i,x+j — X?',x) • (4-2) 



The simplest vortex-related quantity is now [22 



(O tot > = (|n iijX |). (4.3) 
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a) 



b) 



Figure 1: If two vortices pass through one cell, the vortex tracing algorithm must 
decide how to connect them, and this leads to an ambiguity in the length distribution. 



It has the interpretation of the total vortex density, independent of their length. It 
is a perfectly well defined local observable, and any singularities in its finite volume 
behaviour are expected to be related to the thermodynamic singularities of the system. 

However, in the continuum limit (t — f3 — f3 c — > 0) O to t smoothly approaches a 
constant value ~ 0.160, and the interpretation of O to t as a vortex density becomes 
problematic: if we rescale the length scale by the physical correlation length £ ~ t~ v , 
the vortex density diverges as t~ 2u . From the continuum physics point of view, the 
density of vortex lines increases as we increase the resolution at which we look at the 
system, and finally we are measuring only ultraviolet fluctuations. This implies that 
quantities like the number of vortex loops and the length of individual loops are not 
well defined. The same problem arises also in more general continuum theories (see 
25 and references therein). 



These problems might be circumvented if we only monitor sufficiently long loops. 
Thus, we next connect single plaquettes with n^- x ^ to a macroscopic closed vortex 
loop. This steps involves a lot of ambiguity, see Fig. |l|. We will employ two different 
"connectivity" definitions here: 

I. In the maximal definition, in each ambiguous case we make the unique choice so 
that the separate vortex loops are joined together, see Fig. [l].a. 

II. In the stochastic definition, in each ambiguous case the choice is made by throwing 
a dice. This is the procedure most commonly followed in the literature. 

We can thus extract from each lattice configuration a set of "vortex loops" , consisting 
of a number of plaquettes which have been bound together by one of the connectivity 
definitions above. From this set of loops, we measure (using always cubic volumes iV 3 ) 
the following observables: 
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1. "Closed loop order parameter", O c i ose d- For a given configuration, O c i OS ed = 1 if 
there is a vortex which forms a closed non-contractible loop around the lattice with 
periodic boundary conditions (in fact, there have to be at least two such loops, because 
the total topological charge is zero), and O c i osc d = otherwise. We observe that (O c i ose d) 
approaches a step-like function at iV — ► oo, being non-zero in the symmetric phase of 
the XY model. We can also consider a variant of this observable in which we count all 
vortices that extend over a distance > N in any direction (O xyz ) or in some particular 
direction (O x ). 

2. "Volume", O m \- We fit each vortex loop to the smallest possible rectangular 
box ni x n 2 x n 3 . Let O vo \ be the volume of the largest box compared with the total, 
O v oi = nin 2 n 3 /N 3 . We again observe that (O vo i) approaches a step function in the 
infinite- volume limit. 

3. "Mass", O mass . For each configuration, we count the number of plaquettes visited 
by the longest loop L , and normalise by the volume: O mass = (1/-/V 3 ) yj plaq gLo . We 
observe a vanishing result in the broken phase for N — ► oo, a non- vanishing in the 
symmetric phase, but a continuous behaviour at fi c . Thus (O mass ) behaves like the 
magnetization M but with an "inverted" /3-axis. 

4. "Line tension", Tl- If in the idealised case of non-interacting loops, the tension, 
i.e. the mass per unit length, of a vortex is T, then the loop length distribution p(L) 



has the form 12 



p(L) ~ L- 5/2 exp(-LT). (4.4) 

In the interacting case, the line tension Tl is defined by fitting the loop distribution 
p(L) to a similar function, 

p(L)~L- Q exp(-LT L ), (4.5) 

where a and Tl are free parameters. However, in the interacting case Tl should not 
be confused with the tension T given by the free energy per unit length of an isolated 
vortex, and the two do not in general coincide. We observe that the line tension behaves 
similarly to the real (not the inverted) magnetization M. 

For any observable O;, i = {closed, xyz, x, vol, mass}, we define the percolation point 
fii as the point where (Oi) vanishes when fi is increased. Because of the similarities 
in the definitions of the observables, we expect there to be certain inequalities. Even 
in a finite system and in any configuration, the observables satisfy O xyz > O c i ose d, 
O xyz > O x , as well as O vo \ > O mass /3 because each lattice cube has three plaquettes. 
The lattice symmetry also implies that (O x ) > (O xyz )/3, because on the average at 
least every third long loop has to extend into the x-direction. These relations imply 

Pmass ^ Pvoly /^closed ^ ftxyz fix- (^"^v 

It also seems very plausible that in the infinite- volume limit, almost all configurations 
have either O vo \ = or 1 and consequently O xyz > O vo \. We may also conjecture that 
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in the infinite volume limit (O vo \) = (O xyz ), and that (O mass ) > implies (O c iosed) > 
because long volume filling loops have a finite probability to close. Therefore, we may 
expect that 

Z^mass — /^closed — ftvol Pxyz Pxy (^•'O 

and we cannot rule out the possibility that all these points are equal. However, in our 



actual analysis we only rely on Eq. (|4.6|), and use Eq. (4.7) just as a guidance 



5 Results 

We have carried out simulations using both Swendsen-Wang and Wolff cluster algo- 
rithms. We use volumes in the range 4 3 ...200 3 . Individual runs have up to ~ 10 6 
measurements at smaller volumes, ~ 10 4 at the largest ones. 



5.1 Calibration 

As a calibration of our accuracy we first briefly discuss the thermodynamic critical 
properties of the XY model, and show that we can reproduce the known results with 
good precision. 

As an example, we show C tot in Fig. |](a) (cf. Eqs. (|3.3| ), ( fOD ). The critical point, 
defined as the point where C tot obtains its maximum, is fit at relatively small volumes 
4 3 ...48 3 according to Eq. ( |3.6|) , with exponents fixed according to Eq. G3.7T) . We obtain 
f3 c = 0.45430(27), x 2 /dof = 0.19, in agreement with Eq. flSTTp. We can also leave 
both v and u> as free parameters in the fit, but then the errors become quite large: 
v = 0.83(15), oj = 1.16(42), (3 C = 0.45389(62), x 2 /dof = 0.16. However, they are 
consistent with the values in Eqs. (|Q|), ( |3.7| ). As another example, we sit precisely at 
the j3 c given in Eq. ( |3.1| ), and measure the finite volume behaviour of magnetization 
with volumes 20 3 ...200 3 . We find the exponent (3/u = 0.5192(15), in agreement with 
Eq. (|3.7|) . We conclude that our lattice sizes and statistics are sufficient to resolve the 
physical critical properties of the XY model. 



5.2 Volume order parameter 

We now move to our actual topic: percolation related observables. As an illustration 
of a behaviour similar to but not identical with that in Fig. |](a), we first discuss the 
volume order parameter O vo i, determine its percolation critical point j3 vo \, and compare 



it with f3 c in Eq. (|3~T 



In Figs. ^|(b),(c) we show the susceptibility C vo i, as well as its maximum position, 
measured with both types of connectivity definitions. We make an infinite volume fit 
according to Eq. (|3.6| ), with exponents fixed as above, and another one with two free 
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Figure 2: Left: different susceptibilities, Eq. (|3.3| ), as a function of /3; (a) C to t, (b) 
C vo i with the maximal connectivity definition, (c) C vo i with the stochastic connectivity 
definition. Right: the location of the apparent critical point determined from the 
maximum of C. The lines indicate the location of /3 C . 



exponents. With the maximal connectivity definition we obtain 

(3 vol = 0.46184(14), i/ = 0.67, uj = 0.79, xV^of = 2.35, 

Aroi = 0.46074(48), v = 0.66(3), u = 0.12(3), x 2 /dof = 0.64, 

and with the stochastic definition 

Aroi = 0.45565(13), ^ = 0.67, uj = 0.79, xV dof = !- 25 > 



(5.1) 
(5.2) 

(5.3) 
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Figure 3: Left: the reweighted (O c i ose d) as a function of (3. Right: the apparent 
percolation point as obtained from the position of the maximum of C c i ose d (or, for an 
observable only having values 0,1, from the point where (O c i ose d) = 0.5), together with 
the fit in Eq. ( |5.5| ) whereby /3 c i OS ed = 0.45435. The log-log plot shows clearly that the 
power law describes the data well and /Closed > Pc for iV — ► oo. 



(3 ml = 0.45485(37), v = 0.77(8), u = 0.66(19), x7 dof = 0.07. (5.4) 

We observe that the location of /3 vo \ depends on the definition used in constructing 
macroscopic vortex loops in a statistically significant way. The ambiguity inherent 
in this procedure is due to discretization (lattice) effects, cf. Fig. [I]. Moreover, for 
the maximal connectivity definition, the percolation critical point is (3 vo \ ~ 0.461(1), 
far above the thermodynamic critical point, as is already apparent in Fig. 0. Even 
for the stochastic connectivity definition, f3 vol does not coincide with the (3 C , with a 
difference on the level of (2...3)<t depending slightly on the ansatz used for the N — > oo 
extrapolation. For a more precise statement we rely below on Eq. ( |4.6|) and /3 mass . 



5.3 Closed loop order parameter 

Consider then O c i ose d- As conjectured by Eq. (P~7|), /Closed can be even closer to f3 c than 
/3 vo i, and we have therefore carried out a more precise analysis with volumes up to 200 3 . 
The behaviour is shown in Fig. A fit according to Eq. ( |3.6|) to volumes > 16 3 with 
one free exponent gives 

Closed = 0.45435(3), ^ = 0.490(6), x7dof= 9.5/8. (5.5) 

The data is not accurate enough to allow for a good fit with two free exponents. A fit 
with exponents fixed to the thermodynamic values of v and u has x 2 /dof = 53/9, and 
thus a very bad confidence level ~2 x 10~ 8 . 
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Figure 4: The behaviour of (O mass ) at /3 C , compared with the thermodynamic quantity 
M defined in Eq. ( |3.2| ) (more precisely, M/4 in order to have similar magnitudes). 

The central value of /Closed in Eq- fl5.5| ) is now much closer to (3 C than in Eq. (|5.4|). 
However the errors are also an order of magnitude smaller, due to the larger volumes 
used and the smaller range of reweighting in f3. Thus the discrepancy is ~ 6a. The 
disagreement can be very clearly seen in Fig. ^(right). 

We do not consider separately the case of maximal connectivity, since it would in- 
crease /? c iosed moving it further away from (3 C . Likewise, Eq. ( |4.6| ) shows that (3 X and 
(3 xyz cannot coincide with (3 C , either. 

5.4 Mass order parameter 

Since both /3 vol and /3 c i ose d are above /3 C , Eq. (|4.7|) suggests that /3 mass may be very 
close to it. A direct determination of /3 mass is difficult because it turns out that there 
is no clear peak in C mass up to very large volumes, due to strong fluctuations in the 
symmetric phase, and thus it is not easy to obtain finite volume estimates for /3 mass . 
Therefore, we use a different technique to demonstrate that /3 mass > j3 c . 

We stay precisely at /3 C , given in Eq. (|3.1| ), and measure (O mass ) on lattices of sizes 
up to 200 3 . The results are shown in Fig. f|. Recall that we expect a vanishing result 
at /3 mass , and a finite value in the percolated phase (3 < /3 mass . Since (O mass ) behaves 
similarly to magnetization with an inverted /3-axis, we also show the magnetization in 
the same figure. 

A fit to volumes 20 3 . . . 200 3 with a single free exponent gives 0.60/iV a61(4) +0.0393(7), 
with x 2 /dof = 7.3/5. This clearly shows a non-vanishing value in the thermodynamic 
limit. For the magnetization, a fit with a free exponent with volumes 20 3 . . . 200 3 gives 
11 y iV o.5io(9) _ o.002(3), x 2 /dof = 0.49, consistent with 0/u = 0.52 in Eq. Q and a 
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Figure 5: Left: The length distribution of vortices at /3 slightly above /3 C with the 
stochastic definition. Right: the line tension, determined from the fit in Eq. (O) with 
a free parameter a, with both definitions. 



vanishing result at iV — > oo. Thus, while magnetization exhibits the expected critical 
behaviour, from the point of view of O mass we are already in the percolated phase, 
showing that ft mass > ftc- In the light of Eq. fl4.6|), this also confirms that ft vo \ > ft c . 



5.5 Line tension 

Let us finally discuss the line tension. The data is shown in Fig. ^. Since the analysis 
of Tl requires simulations at several values of ft > ft c , we have performed the analysis 
using volumes 48 3 and 64 3 only. We will again monitor ftr L , as well as a critical 
exponent Jt l , defined by 



T L cx \(3-(3 Tl \ 1t ^. (5.6) 

Since this behaviour assumes the infinite volume case, we cannot in practice go too 
close to ft c where the finite volume cuts off the critical fluctuations. At the same time, 
we should be close enough to ft c to be in the critical regime. In practice, we have 
chosen the window ft G (0.46, 0.48) to satisfy these requirements, and have checked the 
volume independence at a couple of ft, see Fig. ||(right). 

Using the stochastic definition, a fit to Eq. ( |4.5|) in this range gives a = 2.11(2), and 
T L as shown in Fig. |](right). Fitting this T L to Eq. (|5.6|) , as well as the one obtained 
with a fixed a = 2.5, we get 



a 



a 



2.11(2), 
2.5, ft TL 



ft Th = 0.4555(2), 
= 0.4558(4), 7Ti 



lTh = 1.62(3), x7dof=0.27, (5.7) 
= 1.67(5), x 2 /dof=2.7. (5.8) 



If ftr L is fixed by hand to a smaller value, 7t l can go up to ~ 1.72. The maximal 
connectivity definition gives clearly contradicting numbers, see Fig. |5|(right), with an 
exponent = 1.80(3). 
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In conclusion, the percolation point (3t l again differs significantly from the critical 
point in Eq. (|3.1|). Moreover, the values we find for 7^ L do not coincide with any of 
the thermodynamic exponents.^ 



6 Conclusions 

It is quite obvious that in discrete spin models even naive geometrically defined perco- 
lation observables show structure close to the thermodynamic phase transition. Indeed, 
the ordered and disordered phases of the spin models are quite well understood, and 
either almost frozen (in the ordered phase), or with almost random fluctuations (in 
the disordered phase), which guarantees that percolation seems to take place. The 
question we have addressed here is whether this rough agreement can be promoted to 
a precise level in the vicinity of the thermodynamic transition point, where the spin 
models show universal critical behaviour, and allow to define a non-trivial continuum 
field theory. 

We have found that in the three-dimensional XY model, the geometrically defined 
percolation transition of the vortex loops is extremely close to the thermodynamic 
phase transition, but does not coincide with it with any of the observables we have used. 
There is a clear discrepancy between the "maximal" and "stochastic" definitions used 
for constructing macroscopic vortex loops. Even with the stochastic definition, a finite- 
size scaling analysis at the thermodynamic critical point, whose location is known very 
accurately from previous studies, shows the absence of any true critical behaviour. We 
of course cannot exclude that it would be possible to bring the percolation transition 
even closer to the thermodynamic one by very carefully fine-tuning the percolation 
criteria; however, this kind of a procedure would have no predictive power. 

We believe that these conclusions apply also to more realistic theories, such as the 
Ginzburg-Landau theory of superconductivity and non-Abelian gauge theories. These 
theories are typically much more complicated than the XY model used here, and there- 
fore it is in most cases very difficult to achieve the level of accuracy needed to confirm 
this in a Monte Carlo simulation. However, we may note that the observation of a per- 
colation transition related to so called Z-vortices in crossover regime of the electroweak 
theory [O] bears a similarity with our conclusion. 

Finally, we would like to emphasise that although geometrically defined percolation 
observables may not reflect the properties of the true phase transition in a rigorous 
sense, the description of the system in terms of the vortex loop degrees of freedom can 
still be useful, if a good observable is chosen. This is certainly true for the tension T 
as discussed in the Introduction, and might also be the case in the non-equilibrium 



6 Our result for the exponent 77^ does not agree with 7t l = 1.45(5) measured in [1J], nor with 
-1t l = 1-50(1) measured in [jl3|. However, the fitted fi ranges there were about 10 Q or more |H| 
times wider than in our case, so that corrections to asymptotic scaling could play a role. 
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dynamics well after a cosmological phase transition, when the vortex loops cease to be 
merely thermal fluctuations and become macroscopic classical objects . 
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